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Introduction 


This  is  the  report  of  activities  carried  out  under 
Contract  F49620-77-C-0031  for  the  calendar  year  1977.  The 
purpose  of  this  contract  is  to  demonstrate  the  feasibility 
of  a new  technique  to  study  unsteady  turbulent  boundary 
layers  experimentally.  The  difficulties  of  making  careful 
measurements  on  a moving  airfoil  or  generating  an  airstream 
that  is  unsteady  and  uniform  are  well  known.  The  former  leads 
to  problems  of  driving  the  airfoil,  of  extracting  the  data, 
as  well  as  correcting  for  transducer  reading  errors  caused  by 
inertial  loads  which  appear  as  spurious  unsteady  forces  or 
pressures.  The  latter  presents  its  own  problems  of  power 
consumption,  complicated  geometries  and  installation. 

The  experimental  studies  performed  under  this  contract 
demonstrate  feasibility  of  a different  type  of  experimental 
technique.  The  principle  underlying  the  technique  is  simple. 

The  circulation  about  an  airfoil  is  controlled  by  the  angle 
of  attack  at  its  trailing  edge.  This  angle  can  be  changed  by 
distorting  the  shape  of  airfoil  or  by  aerodynamic  interference. 
This  interference  is  the  consequence  of  a suitably  shaped 
external  body  located  to  change  the  angle  of  attack  or  the 
effective  airfoil  chamber.  If  the  external  body  moves  nearer 
and  further  from  the  trailing  edge,  the  circulation  varies  in 
a known  fashion.  In  this  way  one  can  control  the  pressure 
gradient  over  the  airfoil.  If  the  external  body  moves  quickly, 
the  pressure  distribution  is  unsteady,  as  is  the  boundary  layer. 
However,  the  airfoil  and  the  measuring  surface  is  and  remains 
stationary.  Thus  with  proper  instrumentation,  pressure  and 
boundary  layer  measurements  can  be  made  with  little  more 
difficulty  than  in  steady  flow. 

Since  the  technique  of  controlling  the  trailing  edge  angle 
of  attack  had  not,  as  far  as  the  authors  know,  been  used 
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previously  in  the  present  form  , the  primary  activities  for 
the  first  year  consisted  of  experiments  designed  to  support 
the  validity  of  the  procedure  and  to  partially  define  the 
characteristics  of  the  boundary  layer  in  the  wind  tunnel  where 
the  experiment  is  to  be  conducted.  For  the  initial  effort/ 
the  steady  boundary  layer  properties  in  a zero  pressure  gradient 
were  to  be  measured.  Ultimately  the  boundary  layer  is  to  be 
developed  on  the  upper  surface  of  a fixed  airfoil,  an  NACA  0012, 
to  be  specific.  Because  it  is  necessary  to  measure  the  pressure 
distribution  over  the  airfoil,  including  near  the  trailing  edge, 
there  is  a possibility  that  some  approximate  data  on  the  un- 
steady Kutta  condition  will  also  be  generated. 

Ultimately  these  studies  could  provide  detailed  informa- 
tion about  unsteady  boundary  layers.  This  information  is  of 
technical  interest  for  three  reasons.  First,  the  unsteady 
characteristics  of  turbulent  boundary  layers  do  not  seem  to  be 
well  understood.  Second,  there  are  numerous  technical  applica- 
tions for  which  the  information  developed  could  lead  to  improved 
device  performance.  Third,  data  of  this  kind  could  help  validate 
results  found  by  computation. 

General  Remarks 

Even  a disinterested  observer  would  be  able  to  note  an 
increased  level  of  activity  in  the  study  of  unsteady  aerodynamics. 
The  most  obvious  manifestation  of  this  interest  is  in  both  the 
AFOSR  meeting  at  the  University  of  Arizona  (1)  and  the  AGARD  Fluid 
Dynamics  Panel  meeting  on  Unsteady  Aerodynamics  (2).  This  in- 
terest may  have  been  stimulated  because  unsteady  aerodynamics  is 
one  of  the  physical  processes  that  govern  or  may  even  limit  the 
performance  of  helicopters,  gas  turbine  engines  and  wind  turbines 
used  as  a power  source. 

_ 

It  should  be  noted  that  A.M.O.  Smith  in  his  Wright  Brothers 
Lecture  "High-Lift  Aerodynamics"  (reported  in  Journal  of 
Aircraft,  1975)  used  this  technique  theoretically  to  study 
high-lift,  multiple  element  airfoils. 
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Another  reason  for  the  increased  interest  may  stem  from 
the  realization  that  many  aspects  of  unsteady  aerodynamics 
may  now  be  studied  through  uses  of  large  scale  automatic 
computation.  Computational  study  is  particularly  attractive 
because  the  difficulties  which  limit  analytical  or  classical 
experimental  studies  do  not  limit  studies  carried  out  by  the 
computer.  The  experimental  data  is  useful  in  supporting 
computational  studies.  The  synergysm  that  results  from  the 
use  of  results  prom  theory,  experiment  and  the  computer  cannot 
fail  to  lead  to  a fuller  understanding  of  a particular  aero-  » 

dynamic  phenomenon. 

Current  Accomplishments 

During  the  last  year  a number  of  preliminary  studies 
were  made.  These  studies,  which  constituted  necessary  pre- 
requisites to  the  main  study,  are: 

1.  Measurement  of  the  characteristics  of  a turbulent 
boundary  layer  on  a zero  pressure  gradient  flat  plate. 

2.  Verification  of  the  ability  to  measure  pressure 
distribution  on  an  0012  airfoil  in  uniform  flow. 

3.  Verification  of  the  ability  to  vary  circulation 
using  the  aerodynamic  interference  near  the  trailing  edge 
of  the  airfoil. 

The  results  from  each  of  these  studies,  which  are  abstracted 
from  Reference  3,  will  be  discussed  in  more  detail  below. 


Wind  Tunnel 


Experimental  Equipment 


All  tests  reported  here  were  conducted  in  the  Wright 
Brothers  Wind  Tunnel  (WBWT)  at  M.I.T.  It  is  a closed  return 
facility  with  a 10  foot  (3.05  meter)  wide,  7.15  foot  (2.29 
meter)  high  elliptical  test  section  15  feet  (4.57  met_r) 
long.  A 2000  HP  motor  provides  wind  speed  up  to  140  miles 
per  hour  (62.5  meters  per  second)  at  atmospheric  pressure. 

Bo undary  Layer  Experiment 

The  two-dimensional  flat  plate  provided  the  surface  on 
which  the  boundary  layer  was  developed.  The  dimensions  of 
the  plate  were  36  inches  (0.9144  meter)  long,  1 inch  (0.0254 
meter) thick  and  20  inches  (0.5080  meter)  wide  at  the  leading 
edge.  The  width  of  the  gap  at  its  trailing  edge  was  slightly 
larger  to  compensate  for  the  streamline  displacement  due  to 
the  boundary  layers  on  the  end  plates.  The  use  of  end  plates 
helped  insure  the  two-dimensionality  of  the  flow  over  the 
flat  plate.  The  dimensions  of  the  end  plates  were  60  x 75  x 
1 in.  (1.5240  x 1.9050  x 0.0254  meter).  The  end  plates  had 
three  static  pressure  taps  and  the  plate  itself  had  five 
static  pressure  taps.  The  leading  edge  of  the  plate  started 
10  inches  (.254  m)  downstream  of  the  48  inch  (121.9  m)  strips 
of  plywood  that  made  up  the  side  walls.  The  leading  edge  of 
the  strip  was  beveled  from  the  outside  at  15°  with  a radius  of 
roughly  1/8  inch  (0.0032  m) . The  bottom  of  the  plate  leading 
edge  was  roughly  a tangent  ogive  of  radius  50.5  inches  (1.282  m) . 
The  tangent  point  was  10  inches  (2.54  m)  from  the  leading  edge. 
The  leading  edge  was  sharp.  The  upper  surface  of  the  plate 
was  Blanchard  ground  and  was  smooth  to  within  35  iricroinches 
(8.98  micrometers).  Static  pressure  taps  were  located  along 
the  center  line  of  the  plate  and  off  the  center  line,  as  shown 
in  Figure  2.  The  end  plates  were  angulated  in  yaw  until  the 
pressure  distribution  was  essentially  flat.  The  resulting 
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gups  were  sealed  with  plasticine.  Figure  2 shows  the  velocity 
distribution  along  the  top  of  the  plate,  normalized  with  re- 
spect to  the  wind  tunnel  "speed  holes".  The  two  cases  shown 
correspond  to  the  final  pitch  angle  setting,  with  and  without 
the  boundary  layer  probe  drive  mechanism.  The  axial  velocity 
ratio  was  relatively  constant  over  the  last  8 inches  (.2  m) 
of  the  plate  at  1.060  ± .002. 

The  velocity  distribution  in  the  boundary  layer  on  the 
flat  plate  was  measured  using  a Flow  Corporation  series 
900  hotwire  anemometer  operated  in  the  constant  temperature 
mode.  The  active  element  of  the  system  was  a platinum  wire 
1 in.  (2.54  cm)  long  and  0.0005  inches  (0.01272  nun)  in 
diameter.  The  hotwire  was  held  by  a cylindrical  probe 
attached  to  a streamlined  arm.  A probe  drive  mechanism  with 
two  degrees  of  freedom  was  attached  to  the  plate.  Vertical 
displacement  of  the  hotwire,  as  indicated  by  a linear 
potentiometer,  was  measured  by  a digital  voltmeter  and  was 
recorded  on  a Hewlett-Packard  model  3400A  two-dimensional 
plotter.  The  hotwire  anemometer  had  been  calibrated  by 
using  the  relationship  between  the  output  voltage  of  the 
linearizer  and  the  velocity  of  free  stream  measured  by  a 
pitot  tube.  The  calibration  constant  was  found  for  the 
probe  positioned  outside  of  the  boundary  layer.  The  probe 
drive  mechanism  and  calibration  constant  were  used  to  de- 
termine the  local  velocity  u (y) . A standard  pitot-static 
probe  was  also  used  to  measure  local  velocities  in  the 
boundary  layer.  When  corrected  using  the  Young  and  Mass 
displacement  effect,  the  two  methods  give  the  same  velocity 
profile. 
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Typical  veloci ty  prof  i los  are  shown  in  Figures  2 and  3 
in  the  form  of  the  inner  variable.  For  these  curves  u 

is  based  upon  Schlicting's  results  (Cf  = .0034)  at 

I r 

Re^  ^ 6,500  (4).  This  leads  to  an  inner  law  of  the  form 

yu 

— » 5.75  1 og . _ + 5.5 

u J 1 0 v 

T 

' 

in  good  agreement  with  the  constants  given  in  Duncan, 

Thom  and  Young  (5).  By  use  of  the  suggestion  of  Clauser  (6) 
who  writes;  using  the  free  stream  velocity,  v^ 


one  can  iterate  graphically  to  obtain  a consistent  fit  of 
A,  B and  Cj.  In  this  way  we  obtained  = .003,  A = 5.75 

and  B = 5.1. 

These  results  are  in  reasonable  agreement  with  those 
reported  by  Fernholz  (7).  Those  results  suggest  the 
turbulent  boundary  layers  generated  in  the  Wright  Brothers 
Wind  Tunnel  will  be  adequate  for  the  purposes  of  the  ex- 
perimental work  planned. 

Ai r f o i 1 

The  0012  airfoil  was  based  upon  coordinates  given  by 
Abott  and  von  Doenhoff  (8).  The  airfoil  was  machined  in 
two  halves,  with  19  pressure  taps  on  each  surface.  The 
pressure  taps  were  connected  to  a Setra  pressure  transducer 
through  a Scannivalve  pressure  scanning  value.  The  wing 
and  its  end  plates  are  shown  in  Figure  4.  Figure  5 shows 
t lie  0012  airfoil  develops  local  velocities  within  about  ± 1% 

I of  those  given  by  Abott.  Hence,  the  airfoil  is  assumed  to 

be  adequate  for  these  tests. 
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Demonstration  of  Principle 
of  Aerodynamic  fill  crfci  tmce 


For  purposes  of  analysis  the  usu<’ 1 assumptions  of 
linearized  aerodynamics  of  an  incompressible  fluid  were 
adopted.  The  local  tangential  perturbation  velocity  can 
be  computed  from  the  integral  given  in  Thwaites  (9) . The 
origin  is  located  at  the  nose  of  the  airfoil. 


o 


here  + 


x 

c 


u 

V 

CO 


e 

a (x) 


is  the  upper  surface 
is  the  lower  surface 
station  of  interest 
airfoil  chord 
perturbation  velocity 
free  stream  velocity 

running  variable  of  integration 

distribution  of  angle  of  attack  or 
camber  over  the  airfoil,  measured 
in  radians 


A circular  cylinder  of  radius  a with  center  located  at 
x = l and  z — — z1  (Figure  6a)  causes  an  angle  of  attack 
distribution  along  the  chord  line  of  the  airfoil. 


a(x)  = tan 


-1 


2 

jjZ-x)2+z,2J  -a2£(*-x)2- z,2J 


2a"  (C-x'z* 
2 


For  each  distance  from  the  trailing  edge  to  the  center  of  the 
cylinder  the  maximum  induced  angle  of  attack  at  the  trailing 
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edge  can  be  computed.  The  locus  position  of  the  center 
of  the  circle  is  shown  in  Figure  6b.  It  is  convenient  to 
normalize  all  lengths  with  respect  to  the  chord;  i.e., 


If  a < 8°  this  expression  can  be  simplified  using  the 
linearizing  approximation 


. -1 
tan  y ± 


3 

3 


= y 


The  term  (~)  in  the  denominator  is 

• z * 

case.  For  simplicity,  we  denote  — 


neglibly  small 

— £ — 
= h and  — = £ , 
c 


for  this 
finally: 


2 ai  (t  - |i  h 

- V+  w 

c 

This  expression  for  the  angle  of  attack  distribution  is  then 

•k 

substituted  into  the  integral.  It  is  evaluated  most  con- 
veniently through  the  variable  change 


* 

For  the  purpose  of  this  calculation,  we  are  neglecting  the 
effect  of  the  airfoil  in  flow  over  the  cylinder.  The  error 
so  incurred  is  not  large  if  the  cylinder  is  not  too  close 
to  the  airfoil. 


X = (1  + cos  0) 


and  the  assumption  that  on  the  airfoil 


ot(0)  = B + B.  cos  0 + ...  + B cos  n0  + 
o l n 


from  which  it  may  be  deduced  (9) 


u_ 

v 


= ± B 


tan  % + B sin  n0 

o 2 ^ — n 


n=l 


The  term  Bq  represents  the  constant  angle  of  attack  and 
represents  curvature  effects. 


Figure  6c  shows  the  distribution  of  a due  to  the 

0.125"  cylinder  so  located  that  a = <s  „ = 8°  at  the  trailing 

edge.  Note  a is  largest  near  the  trailing  edge,  the 

average  angle  of  attack  is  fairly  small,  but  some  camber 

effect  is  apparent.  The  distribution  of  c^  is  compared 

with  the  measured  value  in  Figure  7.  The  difference 

between  the  measured  c and  the  calculated  c , on  the  upper 

P P 

surface  is  partly  due  to  the  inaccuracy  introduced  by  the 

tan  ~ term,  which  becomes  singular  at  x = 0 , and  partly 

due  to  the  separation  of  the  flow  over  the  back  of  the 
cylinder.  The  effect  of  cylinder  diameter  was  tested  using 
three  cylinders.  The  cylinders  had  a length  of  20  inches 
(0.5080  meter)  and  diameters  of  5 inches  (0.1270  m)  , 2.5 
inches  (0.0635  m) , and  1.25  inches  (0.0318  m) . 
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Having  verified  the  principle  that  this  procedure  is 
valid,  a design  for  a support  mechanism  to  move  the  cylinder 
was  started.  It  turned  out  the  primary  structure  in  the 
wind  tunnel  would  not  support  the  needed  dynamic  loads  that 
arose  at  high  speed  cylinder  motion.  Thus,  an  alternate 
procedure  was  studied. 

One  can  agree  that  the  angle  of  attack  at  the  trailing 
edge  is  related  to  the  curvature  of  the  interfering  body. 
Thus,  if  one  made  an  elliptic  body,  one  could  change  the 
local  angle  of  attack  and  camber  by  rotating  the  ellipse. 
This  is  a simpler  dynamic  problem  by  far,  as  far  as  the 
structure  is  concerned.  The  rotating  ellipse  introduces 
a complication  that  the  spinning  may  introduce  circulation 
about  the  axis  of  rotation.  However,  this  additional  angle 
can  be  measured.  Fortunately,  the  potential  flow  past  a 
circle  can  be  transformed  into  flow  past  an  ellipse  through 
a conformal  mapping.  The  ellipse  was  designed  using  the 
confoima 1 transformation 


The  circle  is  in  the  z plane,  its  radius  r = £ = 2 inches 
(0.0508  m)  . We  chose  k =*  1.25  inches  (0.0318  m)  to  provide 
an  ellipse  whose  semi-major  axis  simulated  the  5 inch 
cylinder  and  whose  semi-minor  axis  simulated  the  1.25  inch 
cylinder.  In  the  spanwise  direction  the  ellipse  had 
20  inches  (0.509  m)  length.  The  ellipse  could  be  rotated 
about  t lie  axis,  while  rotating  it  simulated  the  smallest 
and  the  largest  cylinders. 

The  ellipse  was  fixed  to  the  end  plates  below  and  behind 
the  airfoil.  The  details  of  the  several  configurations  for 
which  measurements  were  made  are  given  in  Appendix  A. 
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The  pressure  data  was  converted  into  pressure 
coefficients  for  the  upper  and  lower  halves  of  the  wing. 
The  pressure  transducer  calibration  relations  were  used 
to  find  c^;  i.e. 


c 


P 


loc . 
upper 


output  upper  volts  x 6.7 


in . alcohol 
volts 


Ah  in.  alcohol 


c 

P 


loc. 

upper 


/(output  lower  volts)  - A 6.7 


alcohol 


volts 


Ah  in.  alcohol 


Those  plots  were  made  for  the  airfoil  aligned  with  the  free 
stream  and  for  different  wing-cylinder  combinations  which 
changed  the  angles  of  attack.  On  most  of  the  plots  the 
distribution  of  the  local  coefficients  of  pressure  was 
plotted  versus  percent  of  chord.  Using  the  following 
definition  of  c 
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Figures  8-21  show  the  pressure  coefficient  on  the 
upper  and  lower  surface  of  the  airfoil.  Figure  8 shows 
the  airfoil  alone;  inserting  a quarter  chord  cylinder  has 
the  effect  shown  in  Figure  9.  Here  the  pressure  distribu- 
tion is  measured  when  the  trailing  edge  angle  is  about  8-1/2 
degrees.  Figure  11  is  for  a value  of  2.1  degrees  for  that 
parameter.  Note  the  pressure  distribution  is  relatively 
uniform  in  each  case.  Figures  12-14  show  the  same  data  for 
a cylinder  of  1/8  chord.  Figures  15-17  show  the  same  re- 
sult with  a 1/16  diameter  cylinder  at  the  trailing  edge. 
Since  the  trailing  edge  angle  varies  as  the  square  of  the 
normalized  cylinder  diameter,  a direct  comparison  between 
Figures  11  and  12  show  the  angle  of  attack  at  the  trailing 
edge  dominates  the  flow  process. 

Finally  a series  of  tests  wore  made  with  the  ellipse. 
Figures  18-2.1  show  the  effect  of  the  angle  of  attack  of  the 
semi-major  axis  on  the  airfoil  pressure  coefficient  distri- 
bution . 


Conclusions  and  Recommendations 


The  results  indicate  the  use  of  aerodynamic  interference 
produces  a predictable  effect  on  the  airfoil.  This  effect 
seems  adequate  for  the  purpose  at  hand.  The  turbulent  boundary 
layer  developed  on  a flat  plate  is  consistent  with  that  re- 
ported by  others. 

It  is  recommended  that  boundary  layer  characteristics  be 
measured  on  the  airfoil  in  steady  flow  and  that  preliminary 
measurements  in  the  unsteady  case  be  undertaken. 
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meters 


FIGURE  2 BOUNDARY  LAYER  OVER  FLAT  PLATE 


FIGURE  3 BOUNDARY  LAYER  OVER  FLAT  PLATE 


FIGURE  5 VELOCITY  DISTRIBUTION  OVER  THE  AIRFOIL 
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Figure  6C  a(X)  Induced  by  Cylinder  vs  X 
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Appendix  A 


Schedule  of  Test  Runs 


The  following  is  a list  of  the  data  measured  during  the 
test  program.  Several  runs  were  eliminated  when  various 
errors  in  the  measurements  were  discovered.  For  convenience, 
the  wing  data  has  been  listed  together  with  cylinder  combina- 
tion and  wing  ellipse  combination. 


Wing  Cylinder  Data 


Run 

Number 

l 

h 

d 

a ellipse 

1 

— 

— 

— 

- 

2 

- 

- 

- 

- 

3 

— 

— 

— 

- 

4 

1.175 

0.275 

0.250 

- 

5 

1.300 

0.500 

0.250 

- 

6 

1.300 

0.500 

0.125 

- 

7 

1.175 

0.275 

0.125 

- 

8 

1.209 

0.325 

0.125 

- 

9 

1.175 

0.275 

0.063 

- 

10 

1.209 

0.325 

0.063 

— 

11 

1.300 

0.500 

0.063 

- 

12 

1.209 

0.325 

0.250 

- 

13 

1.175 

0.275 

- 

0 

14 

1.175 

0.275 

- 

45 

15 

1.175 

0.275 

- 

90 

16 

1.175 

0.275 

- 

-45 

17 

1.175 

0.275 

- 

22.5 

18 

1.175 

0.275 

— 

67.5 

25 

1.300 

0.500 

— 

0 

26 

1.300 

0.500 

- 

22.5 

27 

1.300 

0.500 

- 

45 

28 

1.300 

0.500 

- 

-45 

29 

1.300 

0.500 

- 

67.5 

30 

1.300 

0.500 

- 

90 
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